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IceCube is a large neutrino telescope of the next generation to be constructed in the Antarctic Ice Sheet near 
the South Pole. We present the conceptual design and the sensitivity of the IceCube detector to predicted fluxes 
of neutrinos, both atmospheric and extra-terrestrial. A complete simulation of the detector design has been used 
to study the detector's capability to search for neutrinos from sources such as active galaxies, and gamma-ray 
bursts. 
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1. INTRODUCTION 

The successful deployment and operation of the 
AMANDA detector have shown that the Antarc- 
tic Ice sheet is an ideal medium and location for a 
large neutrino telescope. The detection of atmo- 
spheric neutrinos in agreement with expectations 
[ §, established AMANDA as a neutrino tele- 
scope. Searches for neutrinos from Supernova [ 
|j| , dark matter [ |(| , point sources of muon neu- 
trinos [0] and diffuse sources of high energy elec- 
tron [ and muon neutrinos [ have demon- 
strated the physics potential of a deep ice neu- 
trino detector. However, a much larger detec- 
tor is needed to reach a sensitivity required for 
the detection of many predicted neutrino fluxes. 
Current and proposed detectors at sites in the 
northern hemisphere are discussed in [ . Ice- 
Cube is a projected under ice neutrino detector 
consisting of 4800 PMT on 80 strings distributed 
over an area of 1 km 2 (figure [l]) and instrumented 
at a depth between 1400 m and 2400 m. Detailed 
documentation can be found in [ ^|. A surface 
airshower detector consisting of 160 stations over 
1 km 2 augments the deep ice component by pro- 
viding a tool for calibration, background rejection 
and air shower physics. The Monte-Carlo simula- 
tion tools used to determine and optimize the per- 
formance of IceCube are verified with AMANDA 
data. They correctly describe the cosmic ray 
muon flux and the atmospheric neutrino flux [ 

One of the principal objectives of IceCube is 
the detection of sources of high energy neutrinos 
of astrophysical origin. IceCube is sensitive to 
all neutrino flavors over a wide range of energies. 
Muons can be observed from about 10 11 eT^ to 
10 18 eV and beyond. Cascades, generated by v e , 
v e , v T , and v T can be observed and reconstructed 
at energies above 10 13 eV '. Tau events can be iden- 
tified above energies of about a PeV. Interaction 
with the Earth will modulate the neutrino fluxes 
emerging at an underground detector. 

The potential backgrounds of large undcr- 
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Figure 1. Schematic diagram of the arrange- 
ment of the strings of the IceCube detector at 
the South Pole station. The existing AMANDA 
detector and the SPASE air shower array will be 
embedded in the IceCube array. 

ground detectors are downgoing cosmic ray 
muons, atmospheric neutrinos, and the dark noise 
signals detected in photomultiplicrs. The simu- 
lated trigger rate of downgoing cosmic ray muons 
in IceCube is 1700 Hz while the rate of atmo- 
spheric neutrinos (v^ and v^) at trigger level is 
about 300/day. Depending on the type of sig- 
nal to be searched for this background is re- 
jected using direction, energy, and neutrino fla- 
vor. At energies below wlPeV, neutrino astron- 
omy must focus on upward going neutrinos. At 
energies above 1 PeV, the cosmic ray muon back- 
ground disappears while the low energy cosmic 
ray background can be rejected using an energy 
rejection cut. The background contribution from 
dark noise is not significant because the total dark 
noise rate of the optical sensor in situ is expected 
to be less than 0.5 kHz. 

The highest sensitivity for astrophysical point 
sources can be achieved with muons. The muon 
channel stands out for two reasons. 1.) Muons 
allow a very good angular resolution of 0.7° over 
a wide energy range. 2.) The effective volume for 
muons exceeds the geometric volume of the detec- 
tor by factors of 10 to more than 50 depending 
on energy. Due to the long range of high energy 



3 



mucins the interaction of the can take place 
at distances of tens of kilometers outside of the 
detector. We will focus on the muon detection 
because it provides the benchmark sensitivity for 
some of the fundamental goals of high energy neu- 
trino astronomy. 

2. DETECTOR DESIGN 

The detector consists of 4800 photomultipli- 
ers arranged on 80 strings at depths of 1400 to 
2400 m. The strings are arranged in a regular 
spaced grid covering a surface area of 1 km 2 as 
shown in figure [j]. A surface component of 160 
ice Cherenkov tanks provides a detector for air 
showers with an energy threshold of about 1 PeV. 
Each string consists of a 60 optical modules (OM) 
spaced at 17 m. The geometric arrangement and 
the total number of OM and strings is a result 
of MonteCarlo simulations with a variety of dif- 
ferent geometries. The photomultiplier signals 
will be digitized inside the pressure housing [ 
[Tl| . The digitized signals are given a global time 
stamp with a precision of <5ns and transmitted 
to the surface. The digital messages are sent to a 
string processer, a global event trigger and event 
builder. All time calibrations will be automated. 
The geometry of the detector will be known to 
a precision of better than 2 m initially after the 
deployment. It will be calibrated more precisely 
(< 1 m) with light flashers on board the OM 
and with cosmic ray muons. Both methods have 
been applied successfully in AMANDA. High en- 
ergy signals and complex events can be calibrated 
with powerful lasers deployed with the detector. 
The absolute orientation can be calibrated with 
coincidences with the surface air shower compo- 
nent, IceTop, and over longer time scales, with 
the observations of the shadow of the moon. Once 
events are built in the surface DAQ, data will be 
processed and filtered. A reduced data set will be 
sent to the Northern hemisphere for further pro- 
cessing and data analysis on a daily basis. Twenty 
four hour satellite connectivity will be available 
for important messages. Construction of the de- 
tector is expected to commence in the Austral 
summer of 2004/2005 and continue for 6 years. 
The growing detector will be in operation dur- 



ing construction, with each string coming online 
within days after deployment. 

3. FUNDAMENTAL PERFORMANCE 
PARAMETERS 

In detailed simulations we compared the re- 
sponse of the detector to cosmic ray muons, to 
atmospheric neutrinos and to a hypothetical hard 
neutrino spectrum as generated by a shock accel- 
eration mechanism [ [T^| . The event rates gener- 
ated by all three fluxes, normalized to one year of 
on-time, are listed in table [l]. The event rates are 
given at trigger level (minimum of 5 OM signals 
within a section of a string) and for full event re- 
construction with cuts applied for the rejection of 
the cosmic ray muon background, which we refer 
to as "level 2". Unless specified otherwise, basic 
background rejection cuts (level 2) are applied for 
all results shown. With an assumed flux strength 
of El x dN v jdE v = lO^s-^m-V-^eV for 
cosmic neutrinos, a factor of 10 below the current 
best limit to high energy diffuse neutrinos [ 
we expect more than 1000 signal events. At this 
stage, both, the background from atmospheric 
neutrinos and the background from cosmic ray 
muons amount to roughly 10 5 events atmospheric 
neutrino induced muon events is based on [ |lj ) , 
including about 5000 events attributed to prompt 
decays of charmed mesons [ [18). The application 
of an energy cut will quickly suppress remaining 
cosmic ray muon background. 





Trigger 


Level 2 


Cosmic v 


3.3x10 s 


1.1 xlO 3 


Atm v 


8.2 xlO 5 


9.6 xlO 4 


Atmosph. fjj 


4.1 xlO 10 


10 xlO 4 



Table 1 

Event rates are given for 1 year of signal 
and background. The signal expectation cor- 
responds to a source flux of E^ x dN v /dE u = 
10 -7 s _1 cm _2 sr _1 GeV. The expectation for at- 
mospheric neutrino induced muon events is based 
on [ ^lj) and it includes the prompt component 
according to [ [l8||(rqprn). 
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Figure 2. Distribution of the reconstructed zenith 
angle a) Signal from a i?~ 2 -source at trigger level, 
b) with background rejection (level 2, shaded), 
and c) atmospheric neutrino background at level 
2. Event numbers are normalized to one year. 



Figure || shows reconstructed zenith angle dis- 
tributions of atmospheric neutrinos and of the in- 
jected cosmic neutrino flux. Downgoing atmo- 
spheric neutrinos are highly suppressed at the 
level because low energy cosmic ray muons need 
to be rejected with energy cuts. The zenith angle 
distribution of the astrophysical neutrino spec- 
trum is modified little below the horizon. Above 
the horizon the losses are much smaller than for 
atmospheric neutrinos because of the harder spec- 
trum. However, no energy cut has yet been ap- 
plied yet to upgoing muons. 

3.1. Effective Detector Area 

The detector sensitivity can be expressed in 
several quantities. One possible choice is the "ef- 
fective detector area" . It is defined as 



A eff (£ M ,6 M ) = 



detected 



x A„ 



(1) 



where iVg enerat ed is the number of muons in the 
test sample that have an energy at any point 
within the fiducial volume and an incident zenith 
angle 9^. We take the point of closest approach 
to the detector center (which might also lie out- 



side the geometrical detector volume) for the en- 
ergy dependency. Adetected is the number of such 
events that are triggered or pass the cut level un- 
der consideration. 

The effective trigger area reaches one square 
kilometer at an energy of a few hundred GeV. 
Roughly 50% of all triggered events pass the 
"standard selection" (Level 2) , independent of the 
muon energy. Figure ^ shows the effective ar- 
eas for four separate energy ranges as a function 
of the zenith angle of the incident muon tracks. 
The detector will provide an effective detection 
area of one square kilometer for upward moving 
muons in the TeV range. Above 100 TeV the se- 
lection allows for a detection of downgoing neu- 
trinos, i.e. for an observation of the southern 
hemisphere (cos 8 > 0). In the PeV range the 
effective area for downgoing muons is above 0.6 
km 2 , increasing towards the horizon. This means 
that IceCube can be operated as a full sky ob- 
servatory for PeV to EeV neutrinos. 
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Figure 3. The effective area is shown as a function 
of the zenith angle after applying Level 2 cuts for 
four ranges of energy. 



3.1.1. Angular resolution 

The angular resolution of the detector is an 
important quantity for the search for neutrinos 
from point sources. Using the angular resolution 
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Figure 4. Pointing resolution for muons. Shown 
is the median space angle error of the recon- 
structed direction as a function of the zenith angle 
of the incident track. 



of the detector, background events can be more 
easily eliminated by restricting the search to a 
small angular region about the known direction 
of the object under investigation. The median 
angular resolution (median pointing error) at dif- 
ferent zenith angles is shown in figure ^|. 

The angular resolution is about 0.7° at TeV en- 
ergies and approaches 0.6° near the horizon. We 
expect significant improvement of the pointing 
resolution with further development of the recon- 
struction algorithms, in particular from including 
amplitude or waveform information. 



an E~ 2 spectrum. The harder energy spectrum 
of an astrophysical neutrino flux can be used to 
discriminate the atmospheric neutrino spectrum 
from the signal. We use a very simple and robust 
observable, the number of optical modules which 
have seen at least one photon, called the channel 
multiplicity (A c h), as an energy separation cut. 
Atmospheric neutrino induced events have typical 
channel multiplicities (event sizes) of about 30 to 
60 channels. The assumed high energy signal flux 
dominates the atmospheric neutrino background 
above channel multiplicities of about 200. 
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Figure 5. Energy spectra of detected neutri- 
nos for a hypothetical E~ 2 source before (shaded) 
and after energy cut (solid line) , and atmospheric 
neutrinos before the energy cut (dashed). 



4. SENSITIVITY TO ASTROPHYSICAL 
NEUTRINO FLUXES 

4.1. Diffuse fluxes 

Theoretical models of astrophysical fluxes of 
high energy neutrinos are often linked to the 
known flux of very high energy cosmic rays, which 
are believed to be of extragalactic origin, see for 
example |l5, 16, 17]]. Many models have been 
developed that predict the expected flux of neu- 
trinos from the sum of all active galaxies in the 
universe. First we will consider the sensitivity of 
IceCube to detect a generic diffuse flux following 



The simulated source strength of E 2 x 
dN v /dE v = 1 x 10- 7 s- 1 cm- 2 sr- 1 GeV would 
result in an expectation of 74 signal events in 
one year of operation, compared to 8 back- 
ground events from atmospheric neutrinos, after 
applying an optimal energy cut (A c h > 227). 
The background expectation was calculated us- 
ing the prompt charm prediction of Bugaev et 
al.[ |jgj| , according to which, prompt decays of 
charmed mesons contribute 80% to the final at- 
mospheric sample. After three years of opera- 
tion an overall flux limit of E 2 x dN v jdE v = 
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4.2 x 10~ 9 s _1 cm _2 sr _1 GeVis obtained. This is 
more than two orders of magnitude below the cur- 
rent limit obtained with AMANDA. This sensi- 
tivity would improve by a factor of 2 when using 
the charm prediction of TIG [ |l9) . It should be 
noted that the sensitivity for diffuse v e - and in- 
fluxes may improve the overall sensitivity signifi- 
cantly. The energy spectra of the incident signal 
and background neutrinos are shown in figure |5|. 
The applied energy cut results in an detection 
threshold of about 200 TeV. The sensitivity ob- 
tained after one year of data taking is already well 
below the Waxmann and Bahcall diffuse bound [ 

Apart from the generic case of a E 2 spectrum, 
which is typical for scenarios that involve meson 
production in interactions of Fermi accelerated 
cosmic rays with matter, we give the sensitivity to 
two other examples of diffuse models. Mannheim, 
Prothcroc and Rachcn have calculated an upper 
bound on the diffuse neutrino flux arising from 
photohadronic interactions in unresolved AGN 
jets in the universe. Their flux bound is shown 
in figure | labeled MPR, as a function of energy. 
Also shown is a model by Stecker and Salamon 
[ [TtJ for proton interactions on the UV thermal 
photon field in AGN cores, labeled S&cS. When 
testing the sensitivity of IceCube to these specific 
models we found that IceCube would be sensitive 
to a flux of 2% of the MPR flux, and 2.3 • 10~ 3 
of the Stecker /Salamon flux. Also shown is the 
GRB flux prediction by Waxmann and Bahcall [ 
(dash-dotted line). Here, the sensitivity limit, 
defined as a potential exclusion at 90% C.L., was 
calculated for an observation of 500 bursts (ex- 
pectation for one year) in the northern sky (see 
section [D]) . 

4.2. Sensitivity to Point Sources 

An excess of events from a particular direc- 
tion in the sky indicates the existence of a point 
source. The angular resolution of 0.7° allows us- 
ing a search window of one degree radius that 
greatly reduces the background, while retaining 
a large fraction of the signal. In order to re- 
ject the remaining small number of atmospheric 
neutrinos we add a soft energy cut by requir- 
ing a channel multiplicity of A c h=30. The ap- 
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Figure 6. Expected sensitivity of the Ice- 
Cube detector. The thick solid line indicates the 
90%c.l. limit setting potential for an E~ 2 type 
spectrum for a time period of three years (see 
text for details). The thin solid line illustrates 
the GRB sensitivity for the reference spectrum [ 
0. 



plied energy cut suppresses atmospheric neutri- 
nos of energies mostly below 1 TeV allowing a es- 
sentially background free detection potential for 
point sources. With these parameters we find 
an average flux upper limit of E 2 x dN v /dE v = 
5.5 x 10~ 9 s _1 cm _2 GeV after one year of data 
taking. After five years of operation the sensi- 
tivity level would reach E 2 x dN v jdE v ~ 1.7 x 
10" 9 s^cm^GeV . 

4.3. Gamma Ray Burst Sensitivity 

The origin of the high-energy gamma-ray burst 
emissions seen by orbiting satellites is still unde- 
termined. If the acceleration mechanisms in these 
objects involve hadronic processes then the ob- 
served fluxes of gamma-rays should be accompa- 
nied by neutrinos. The detection of these neu- 
trinos is made easier by using the time stamp 
and direction provided by the satellite observa- 
tion. This small duration of gamma ray bursts 
allows for an essentially background-free search 
opportunity. Waxman and Bahcall [ [l5| calcu- 
lated the expected flux of neutrinos from the sum 
of all gamma-ray bursts by assuming that the ac- 
celeration mechanism is hadronic and is responsi- 
ble for the observed flux of cosmic rays. The en- 
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ergy spectrum predicted by Waxman and Bahcall 
is shown in figure @. A diffuse search for gamma- 
ray bursts involves summing over the observation 
time and spatial search windows for many sepa- 
rate bursts. 

We considered at the moment an observation 
only over the northern sky, where the search will 
not be limited to higher energies due to downgo- 
ing CR muon background. A downgoing search 
would be similar, only it would require a harder 
energy cut. From the remaining rate of 500 bursts 
per 2tt sr and year we would expect 13 neutrino in- 
duced up-going muons at 0.1 background events. 
For this analysis we used a hypothetical observa- 
tion duration of 10 seconds and a spatial search 
cone of 10° centered about the direction of the 
GRB. The upper limit to the associated GRB flux 
would be 20% of the flux calculated by Waxmann 
and Bahcall. 

5. OTHER SCIENCE OPPORTUNITIES 

5.1. Dark matter 

If Weakly Interacting Massive Particles 
(WIMPs) make up the dark matter of the uni- 
verse, they would also populate the galactic halo 
of our own Galaxy. They would get captured 
by the Earth or the sun where they would an- 
nihilate pairwise, producing high-energy muon 
neutrinos that can be searched for by neutrino 
telescopes. A favorite WIMP candidate is the 
lightest neutralino which arises in the Minimal 
Supersymmetric Model (MSSM) . The typical en- 
ergy of the neutrino-induced muons would be of 
the order of « 25% of the neutralino mass. The 
predicted muon rates from WIMPs annihilating 
in the Sun range up to \Q A km~ 2 per year at 
wlOOGeV and more than a 10 4 fcm -2 up to ener- 
gies of «1 TeV. Current limits lie at fluxes of sev- 
eral thousand events at energies above 100 GeV. 
Simulations indicate that IceCube should reach 
sensitivities below 50 muon events per year for 
WIMPs from the sun [ ^0) . Analysis methods for 
WIMP searches be optimized for low energies. 
IceCube could play a complementary role to fu- 
ture direct detection experiments (like CRESST 
or GENIUS) for annihilation in the Earth, and 
may even have a slight advantage over direct de- 



tection experiments for certain low-mass WIMP 
models and annihilation in the Sun. 

5.2. Cosmic rays and airshowers 

Combined with the 1 km 2 surface detector 
IceCube can do unique coincidence and anti- 
coincidence measurements with high energy air 
showers. In addition to providing a sample of 
events for calibration and for study of air-shower- 
induced backgrounds in IceCube, the surface ar- 
ray will act as a partial veto. All events gener- 
ated by showers with E > 10 15 eV can be ve- 
toed when the shower passes through the surface 
array. In addition, higher energy events, which 
are a potential source of background for neutrino- 
induced cascades, can be vetoed even by showers 
passing a long distance outside the array. The 
IceCube-IceTop coincidence data will cover the 
energy range from below the knee of the cosmic- 
ray spectrum to > 10 18 eV. Each event will con- 
tain a measure of the shower size at the surface 
and a signal from the deep detector produced by 
muons with E > 300 GeV at production. At the 
high elevation of the South Pole, showers will be 
observed near maximum so that measured shower 
size provides a good measure of the total shower 
size. The combined measurement of the muon- 
induced signal in IceCube and the shower size at 
the surface will give a new measure of primary 
composition over three orders of magnitude in en- 
ergy. In particular, if the knee is due to a steep- 
ening of the rigidity spectrum, a steepening of the 
spectrum of protons around 3 • 10 15 eV should be 
followed by a break in the spectrum of iron at 
8 ■ 10 16 eV. The method to measure the composi- 
tion has been developed and applied successfully 
with AMANDA and the surface airshower array 
SPASE-2 [||. 

5.3. SUPERNOVAE 

Although the MeV-level energies of supernova 
neutrinos are far below the AMANDA/IceCube 
trigger threshold, a supernova could be detected 
by observing higher counting rates of individual 
PMTs over a time window of 5-10 s. The en- 
hancement in rate of one PMT will be buried in 
dark noise signals of that PMT. However, sum- 
ming the signals from all PMTs over 10 s, signifi- 
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cant excesses can be observed. With background 
rates more than 10 times lower than ocean exper- 
iments, IceCube has the potential to see a super- 
nova and to generate an alarm signal. 

5.4. Cascades and very high energies 

Due to its large and regularly instrument vol- 
ume, IceCube has a high sensitivity for cascades 
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generated by v e 



and v T . Cascades have 



a superior energy resolution because all energy is 
deposited on a small volume of «10 metres di- 
ameter, and the amount of Cherenkov photons 
scales linearly with the deposited energy. At en- 
ergies above a few tens of TeV the detector be- 
comes fully efficient to cascade detection with 
an effective volume comparable to the geomet- 
ric volume of 1km 3 . At energies above 1 PeV 
chances increase to distinguish tau events from 
the J/ e -induced cascades. Extremely energetic 
muon events reach effective areas well beyond the 
geometric area as indicated in figure 0. Effective 
areas at 10 18 eV are estimated to reach more than 
2.5 km 2 . Reconstruction algorithms for the more 
complex event topologies such as from tau events 
and for extreme energies beyond 10 18 eV are not 
developed yet. However, the high resolution and 
high dynamic range data acquisition system of 
IceCube will provide rich information to extract 
the physics from such events. 
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